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We investigate the phase structure of hot and/or dense QCD using the Schwinger-Dyson 
equation (SDE) with the improved ladder approximation in the Landau gauge. We show that 
^ ' the phase transition from the two-flavor color superconducting (2SC) phase to the quark- 

D ' gluon plasma (QGP) phase is of second order, and that the scaling properties of the Majorana 

pL^ , mass gap and the diquark condensate are consistent with mean field scahng. We examine 

the effect of the antiquark contribution and find that setting the antiquark Majorana mass 
equal to the quark one is a good approximation in the medium density region. We also 
study the effect of the Debye screening mass of the gluon and find that ignoring it causes 
the critical lines to move to the region of higher temperature and higher chemical potential. 



(N 
> 

(N 

^ . §1. Introduction 

CN ' The dynamics of quantum chromo dynamics (QCD) are very rich, and dynam- 



ical chiral symmetry breaking is one of the most important features of QCD. In 



\ hot and/or dense matter, chiral symmetry is expected to be restored (see, e.g. 

Refs. andEj)). Furthermore, in recent years, there has been great interest in the 
Oh. phenomenon called "color superconductivity", which occurs after chiral symmetry 

(— I \ restoration at non-zero chemical potential in the low temperature region (for a re- 

cent review, see, e.g., Ref. i3j). Therefore, exploration of the QCD phase diagram, 
including the color superconducting phase, is an interesting and important subject 
' for the purpose of studying not only the mechanism of dynamical symmetry breaking 

' but also its phenomenological applications in cosmology, the astrophysics of neutron 

stars and the physics of heavy ion collisions .^^'^ 

The phase structure of QCD at non-zero temperature with zero chemical po- 
tential has been extensively studied with lattice simulations, but the simulations at 
finite chemical potential has begun only recently and these simulations still involve 
large errors [see, e.g., Ref . (1| and references cited therein]. Thus, it is important 
to investigate the phase structure of QCD in the finite temperature and/or finite 
chemical potential region by various other approaches. 

In various non-perturbative approaches, the approach based on the Schwinger- 
Dyson Equation (SDE) is one of the most powerful tools [for a review, see, e.g., 
Refs. 5J and J5)]. From the SDE with a suitable running coupling at zero tempera- 
ture and zero chemical potential, the high energy behavior of the mass function is 
shown to be consistent with the result derived from QCD with the operator product 
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expanssion and the renormalization group equation. When we use the SDE at non- 
zero chemical potential, we can distinguish the Majorana mass of the quark from 
that of the antiquark by introducing the on-shell projectors of the free quark and 
the free antiquark, which are useful especially in the high density region, where the 
antiquarks decouple. Furthermore, we can include the effect of the long range force 
mediated by the magnetic mode of the gluon, which may have a substantial effect, 
even in the intermediate chemical potential region, as in the high density region.^''^ 
In the SDE analysis, the phase structure of QCD in the finite temperature and fi- 
nite chemical potential region have been investigated, concentrating on the chiral 
symmetry restoration.'^'GSJElEJEJOJ in a previous work,^ solving the SDE with 
the full momentum dependence included, we studied the phase transition from the 
hadron phase to the two-flavor color superconducting (2SC) phase in the region of 
finite quark chemical potential at zero temperature. It was shown that the phase 
transition is of first order and the existence of the 2SC phase decreases the critical 
chemical potential at which the chiral symmetry is restored. 

In this paper, we extend the previous work to non-zero temperature and solve 
the coupled SDE for the Majorana masses of the quark and antiquark {A'^) 

separately from the SDE for the Dirac mass (B) in the low and intermediate tem- 
perature and chemical potential region. The true vacuum is determined by com- 
paring the values of the effective potential at these solutions. In several analyses 
carried out to this time, the SDE without the antiquark Majorana mass (A'^) has 
been used to estimate the order of the Majorana mass gap at intermediate density 
{fi ~ 300 MeV).l3nSHI3 Here, we investigate this antiquark effect in the intermedi- 
ate density region by comparing the results obtained from the SDEs in the following 
three cases: (case-1) coupled SDEs for the quark and antiquark Majorana masses; 
(case-2) SDE for the quark Majorana mass with the antiquark Majorana mass set to 
zero, which is valid in the high density regionj^'^'^ (case-3) SDE for the Majorana 
mass with the antiquark Majorana mass set equal to the quark one, as in several 
anlayses carried out using models with the contact 4- Fermi interaction [see, e.g., 
Refs.IIHl)iIlI)J2ni)JlIl)J211)JlSl)]- We solve these three types of SDEs for the Majorana 
masses and discuss the importance of the antiquark contribution in the region of 
small chemical potential. We find that the antiquark mass is of the same order as 
the quark mass in the low and medium density region (fi < 600 MeV, where /i is the 
quark chemical potential), and setting Z\+ = A~ is actually a good approximation 
for investigating the phase diagram, the quark Majorana mass gap and the diquark 
condensate. Furthermore, we investigate the effect of the Debye screening mass of 
the gluon by comparing the results from the SDEs in the case of zero Debye mass 
with those in the case of non-zero Debye mass. We find that omission of the Debye 
mass of the gluon causes the critical lines in the phase diagram to move to the region 
of the higher temperature and higher chemical potential. 

This paper is organized as follows. In ^ we summarize the quark propagator, 
the gluon propagator and the running coupling that we use in the present analysis. 
Several approximations of the quark propagator are made. We also give formulas for 
calculating the chiral condensate and the diquark condensate. In ^ we present 
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the effective potential for the quark propagator and then derive the Schwinger- 
Dyson equation as a stationary condition of the effective potential. Section |1] is 
the main part of the paper, where we give the results of the numerical analysis of 
the Schwinger-Dyson equation. Finally, we give a summary and discussion in ^ In 
the appendices we summarize several intricate expressions and formulas. 

§2. Preliminaries 

In this section we present the quark propagator, the gluon propagator and the 
running coupling that we use in the numerical analysis. Using the imaginary time 
formalism, we extend the formalism of Ref. 115(1 to non-zero temperature [for the 
real time formalism, see, e.g., Refs. 1241 ], Below, is related to the Matsubara 
frequency at non-zero temperature T as = (2n + l)7riT. In ^2.11 we introduce the 
eight-component Majorana spinor (Nambu-Gorkov field) and give the general form 
of the full quark propagator as a matrix in the Nambu-Gorkov space. The gluon 
propagator with a screening effect is presented in H2.21 We also give the explicit 
form of the running coupling in the analysis. We give a formula to calculate the 
quark-antiquark condensate and the diquark condensate in ^2.31 

2.1. Nambu-Gorkov field and quark propagator 

In the present analysis, we regard u and d quarks as massless, but we consider 
the current mass of the s-quark to be large enough that the strange quark can 
be ignored in the formation of the diquark condensate; we assume that the color 
superconductivity is realized in the 2-flavor color superconducting (2SC) phase,!^ 
where the color symmetry SU{3)c is broken down to its subgroup SU{2)c.*^ 

Because we are interested in the phase structure of QCD, including the color 
superconducting phase, it is convenient to use the eight-component Majorana spinor 
(Nambu-Gorkov field) instead of the four-component Dirac spinor. The Nambu- 
Gorkov field is expressed as 

where C = ij^'y^. In the Nambu-Gorkov basis, the inverse full quark propagator is 
expressed aa^^ 

■a (^^-i^f iPo + f^h''-l-p-Bip) A{p) 

(2-2) 

B{p)^^ = Br{p)5'^ {5"^ - S^^d^"^) + 53(p)5*-',5"3<5^3 ^ 
MpTi' = e,,e^''j,[A+ip)A; + A-ip)A-] , 

*' We believe that the color flavor locked (CFL) phasd^ is realized in the high density region, 
which is beyond the region of our analysis, and that the 2SC phase is more stable than the CFL 
phase near the chiral restoration point. 
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= -e,,e^'%[A+{p)A^ + zA-(p)yl+] , (2-3) 

where i and j are the flavor indices and a and b are the color indices. We have chosen 
the 3-direction as the direction of color symmetry breaking in such a way that the 
quarks carrying the first two colors make a pair. Here A~ and are the on-shell 
projectors for the quark and antiquark in the massless limit [for discussion of the 
massive case, see e.g., Ref. 177))]: 

Note that in our analysis, the chemical potential fi is introduced with respect to 
the quark number. Now, the full quark propagator includes four scalar functions, 
Bi and B^, corresponding to the Dirac masses responsible for the chiral symmetry 
breaking, and and A~ , corresponding to the Majorana masses responsible for 
the color symmetry breaking. As shown in Refs. 11,1 and 1151 . the Dirac masses Bi 
and i?3 obey the constraint 

Im [B,{p)] = -Im m-p)] , {i = l, 3) (2-5) 

while, as shown in Eef. I15|l . the Majorana masses A'^{p) and A^{p) are real and 
even functions of po^ 

A^ip) = A^{-p) = [A^ip*)Y . (2-6) 

By taking the inverse of the expression in Eq. H2-2() , the full quark propagator is 
obtained as 

-^^ (r.\ = ( ^+'^^ -{{po + fJ^h'-1-p-B(j,)}-'Aip)Rz'(j^) \ 

^^F[P) ^ _^^p^_^^^o_^.^_B^_p-^}-iA{p)R-Hp) Rz\p) J ' 

(2-7) 

where 

R+{p) = {{po + m)7° - 7 • P - B{p)} - A{p){{po - ^)7° - 7 • P " B{-p)}-'A{p) , 
R^{p) = {{po - m)7° - 7 • P - B{-p)} - A{p){{po + ^)70 - 7 • P - B{p)}~^A{p) . 

(2-8) 

2.2. Gluon propagator and running coupling 

In previous analyses of the phase structure carried out using the SDE with the 
improved ladder approximation, e.g., in Refs. !9| and [Tip . a gluon propagator with 
the same form as that at T = /i = was used. However, in a hot and/or dense 
medium, the gluon is subject to medium effects and thereby becomes effectively 
massive. In this paper, we include the Debye screening mass of the electric mode 
of the gluon propagator through the hard thermal/dense loop approximation and 
ignore the Meissner masses, as done in Ref. llSp. However, the form of the magnetic 
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mode used in, e.g., Refs. Ej), IBJ and [T5|) cannot be extrapolated to zero chemical 
potential and zero temperature, because the fc|-term is omitted. In the present 
analysis, therefore, we drop the Landau damping of the magnetic mode and include 
the /cl-term. 

It should be noted that the wave function renormalization for the quark prop- 
agator must be kept equal to 1 to ensure consistency with the bare quark-gluon 
vertex adopted in the improved ladder approximation.!^ At zero temperature and 
zero chemical potential, we can show that the SDE with the ladder approximation 
in the Landau gauge does not result in wave function renormalization (see, e.g., 
Ref. ISI))- In the present analysis, therefore, we take the Landau gauge for the gluon 
propagator, while assuming wave function renormalization for the quark propagator 
to be 1, even at non-zero temperature and/or non-zero chemical potential, as was 
done in Ref . Ej) . 

The explicit form of the gluon propagator that we use in this paper is given by 



(2-9) 



where /i;4 = —ik<^ = lirnT (with integer n) and M^, is the Debye mass of the gluon. 
In the hard thermal/dense loop approximation, the Debye mass is given by^^' 



M|, = as 



-{2N, + Nf}T^ + -^n^ 



TT 



(2-10) 



where Nc = Nj = 3 and the value of as is set to the value at an infrared scale, 
as{Ef) [see Eqs. (|2-13|) - (|2-15|) ]. because we are interested in the region in which 
both the temperature and chemical potential are lower than the infrared scale, i.e., 
T, |U < Ef. The quantities 0^f^l{i = 1,2) are the polarization tensors, defined by 

r>(i) - P-L I • p» (2) _ {u ■ kf 

(u • kf - /t2 ' (u . A:)2 - fc2 ' ^ > 

where 

_ _ k^k,^ _ k^ky ^ k^Uy + u^ky u^u^ ^ 2 ir. 

'^'^"^^''^ A:2 ' ^M^'- ^2 ^.fc ^ (u-kf^ ■ ^ > 

The Lorentz four-vector = (1,0,0,0) in the gluon propagator reflects the explicit 
breaking of Lorentz symmetry due to the existence of a temperature and/or chemical 
potential in the rest frame of the medium. 

In the improved ladder approximation,!^ the high-energy behavior of the run- 
ning coupling is determined by the one-loop renormalization group equation derived 
in QCD, and the low-energy behavior is suitably regularized. As a result, at zero 
temperature and zero chemical potential, the high-energy behavior of the mass func- 
tion derived from the SDE with the running coupling is consistent with that derived 
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from the operator product expansion (OPE).^ In the present analysis, we use the 
following two types of running coupling^'^''^^''^ to check the dependence of the 
results on the methods for the regularization of the infrared behavior: 

(I) -(^) = ^ = II^^-/(M/,^e), (2-13) 



f{t,tf,tc) 



f \ if t>tf, 

tJ + 2tj(tf~t,) 1* tf>t>tc 



tf ' 2t 



(II) a.iE) = ^ = ^ , (2.14) 

^ ^ ' Air llNc-2Nf max(t,t/) ' ^ ^ 



where 



£;2 E'2 

* = 1"7^' */ = ln^' tc = ^^-r^, (2-15) 

^^qcd ^^qcd ^^qcd 

with E being the energy scale, ^qcd the characteristic scale of QCD, and Ej and 
Ec the infrared cutoff scales introduced to regularize the infrared singularity. Note 
that the value of the running coupling in the infrared region needs to be sufficiently 
large in order to involve the dynamical chiral symmetry breaking!^ Here we use the 
running coupling (I) for tf = 0.5 and tc = —2.0, around which the various physical 
quantities for the case T = jj = are very stable with respect to changes in 
The running coupling (I) was used in Refs.lH} and II 1)1 . while the running coupling 
(II) was used in our previous work.i^^ Therefore, we investigate the dependence of 
the result on the regularization methods by adopting several values oft/ for both (I) 
and (II). As discussed in the previous subsection, we assume that the mass of the 
strange quark is large enough that we can ignore the s-quark in the formation of the 
diquark condensate of u and d quarks. On the other hand, it is natural to assume 
that the current mass of the s-quark is smaller than Aqcd- In such a case, the effect 
caused by the s-quark should be included in the running coupling. Thus we set 
Nf = 3 and A'^c = 3 in the running couplings H2-13() and (|2-14|) . 

2.3. Condensates 

In this subsection we give formulas to calculate the chiral condensate and the 
diquark condensate at non-zero temperature. Because we use the imaginary time 
formalism, the formulas for the condensates are obtained by replacing the integra- 
tions over p4 of the expressions at T = in Ref. I15|l with the Matsubara sums as 
/ ^ ~^ X^n • ^^i- I15|l , the logarithmic divergence of the integral is regularized 
by introducing the ultraviolet (UV) cutoff yl4 for the four-dimensional momentum 
as —p^ < A^. In the present analysis, we introduce the UV cutoff A for the spatial 
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momentum as \p\^ 
given by 



< A?. Then, the general expression of the chiral condensate is 



(2^ 



tr[5 



(2-16) 



where 1 17) is the state at non-zero temperature and non-zero chemical potential, and 
the trace is taken in the spinor, flavor and color spaces. Summations over the color 
index a and the flavor index i are implicitly taken on the left-hand side of Eq. ()2-16|) . 
As in the case T = 0, the formula for the chiral condensate of the charge conjugated 
quarks is the same as that of the quarks given above^ (i7|['i/'c]a['0c]i'(O)|l7)yi = 

The general form of the diquark condensate is 



(2vr)^ 



tr[e(^)e(^)5i.i275] 



(2-17) 

where e^'^^ and e*^-^^ are antisymmetric matrices in the color and flavor spaces, respec- 
tively: 



{c)\ab 



(/) 



}ij 



(2-18) 



The diquark condensate of the charge conjugated quarks is same as that of the 
quarks, except for sign.15' {f2\{e'^ eab3MBtCl5[i^c]''{0mA = -me'hab3)[i^'^]'}Cj5 



ijjj{0)\f2) A- The explicit forms of the above condensates written in terms of Bi and 
are given in Appendix EI 

In the improved ladder approximation at zero temperature and zero chemical 
potential, the high-energy behavior of the mass function is consistent with that 
derived from the OPE. The chiral condensate calculated from the mass function was 
shown to obey the renormalization group evolution derived from the OPE (see, e.g., 
Refs. i5j)). Then, we identify the condensates, which are calculated with UV cutoff 
A, with those renormalized at the scale A in QCD. Therefore, we scale them to the 
condensates at 1 GeV, using the leading renormalization group formulas. Note that 
at non-zero temperature, the integral over is converted into a Matsubara sum, and 
the cutoff A is introduced for the spatial momentum, not for the four-dimensional 
momentum, as in Ref. 115(1 . Nevertheless, the relations between the condensates at 
the scale A and those at the scale E have the same form: 

'osiA)' 



(r2|V>f(o)|i?)j 



{f2\{e'^eab3)[^l^^]'^Ci^m^) 



(2-19) 



{n\{e'^eab3)[i^^]1Ci^'j{0mA 



where 



9C2iF) 
llNr - 2N 



f 



llNr 



2Nf 2Nc 



(2-20) 



(2-21) 
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§3. Effective potential and the Schwinger-Dyson equation 

In this section we present tlie effective potential for tlie quark propagator and 
tfien derive tlie Schwinger-Dyson equation (SDE) as a stationary condition of the 
effective potential. 

In the present analysis, we are interested in the differences among the values 
of the effective potential corresponding to the 2SC vacuum, the chiral symmetry 
breaking (xSB) vacuum and the trivial vacuum. Then, as explained in ^2.11 we 
assume that the current mass of the strange quark is large enough that it can be 
ignored in the valence quark sector, and we include only u and d quarks in the 
effective action. In other words, we include the strange quark only quark in 

the present analysis. Then, the effective action for the full quark propagator Sp is 
given bj^^ 

nSp] = \ (^-iTrLn(5^i) - iTTisf'^Sp) - ir2vi[SF]^ , (3-1) 

where TV and Ln are taken in all the spaces and /2Pi['S'f] represents the contributions 
from the two-particle irreducible (with respect to the quark line) diagrams. Here, 

S^p"* is the free quark propagator. The factor 1/2 appears because we use the eight- 
component Nambu-Gorkov spinor. In the high temperature and high density region, 
the one-gluon exchange approximation is valid, because the coupling is weak. In 
the present analysis, we extrapolate this approximation to intermediate temperature 
and intermediate chemical potential, and include only the contribution from the 
one-gluon exchange diagram in F^piiSp]'- 

r2m[SF\ = -\mSf ■ igr^ ■ Sf ■ ign ■ D^u) , (3-2) 
where is the quark-gluon vertex in the Nambu-Gorkov basis defined as 

= ( ^ 0^^ -^'^{T^f ) • ^^-^^ 

From the effective action ()3-l() . the effective potential in the momentum space can 
be written as 

v[Sf] = -r[SF]/ 1 d^x 

\^T. I (^^ det{SF{p)} -tv{sp-\p)SF{p)}^ 



2' 

n 



(2vr)3 Z-fJ (2vr) 



hiiSFip) ■ igr^ ■ SF{q) ■ igm ■ iD^!^{p - q) 

(3-4) 
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where In, det and tr are taken in the spinor, color and flavor spaces. 

The SDE is obtained as the stationary condition of the effective potential 
{5V[Sf]/6Sf = 0): 

Sp^ = SP'' - {igr'^ ■ Sf ■ igm ■ D^^ . (3-5) 

In this paper we investigate the effect of the antiquark contribution by considering 

the following three different types of the SDEs: (case-1) SDE with including the 

quark and antiquark Majorana masses properly; (case-2) SDE with the antiquark 

Majorana mass omitted, as an approximation commonly applied in the high density 

regionP'^'^ (case-3) SDE with the quark Majonara mass and the antiquark one 

set equal, as in analyses carried out using models with the contact 4-Fermi interac- 
^ p) m JH) JH 

case-1 

The SDE l|3-5() leads to the following four coupled equations for the four scalar 
functions Si, B^, and A~: 

B,{p) = tY, j -0^\^asD,M-pMl^T^SFii{q)YT^5^^] , 

(3-6) 
(3-7) 

pMrT''SFi2{q)Y{T^fA-^^e^'h^f^ , 

(3-8) 

p)tr[^^^T^SFi2{q)Y{T^)^A;^^e^'h^f^] , 

(3-9) 
where 

= 5"^ - 5"-^ 6^^ , {4'^}"^ = S^'^d^^ , (3-10) 

and the traces are taken in the spinor, flavor and color spaces. The quantity = 
as{E) on the right-hand sides of Eqs. (|3-6|) - (|3-9|) is the running coupling defined in 
Eqs. (|2-13|) and (|2-14|) . As in R,efs.l!Hl.[TT|l and lT^ . we use the angular averaged form 
E'^ = —p^ — q^ as the argument of the running coupling for T > and /i > in the 
present analysis. The explicit forms of the coupled equations in Eqs. (|3-6j) ~ (|3-9|) are 
given in Appendix ^ The SDE with A~ distinguished from at non-zero is 
analyzed in Refs. IT^ and •> 1 ) . where the analyses are done at zero temperature. 

case-2 
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The SDE (|3-5|) with = leads to three coupled equations for the three scalar 
functions Bi, and . The equation for A~ is given by 

(3-11) 

The equations for Bi and B^ are given by setting = in Eqs. (|3-6|) and (|3-7|) . 
The SDE with = is considered to be valid in the extremely high density region 
and is used in many works at high densityP''^'^ 

case-3 

The SDE H3-5|) with A~ = A~^ leads to three coupled equations for the three 
scalar functions Bi, Bj, and A^{= A'^). The equation for A~ is given by 

(3-12) 

The equations for Bi and B^ are obtained by setting A~ = in Eqs. 1)3 •6() and 
H3-7() . respectively. This type of the Majorana mass gap is used in several analyses 
carried out using models with the local 4- Fermi interactioiP^'C^''^''^''^''^ as well 
as the SDE.EHJ 

As we see in Eq. H2-6() . the Majonara masses A^{p) and A~{p) are real, even 
functions of pq, while in general the Dirac masses Bi{p) and B^[p) are complex 
functions. Equation H2-5|) is the constraint on the imaginary part, while no constraint 
is obtained on the real part from general considerations. However, as shown in 
Refs. 11) and .15) for the case A'^{p) = A~{jp) = 0, the structure of the SDE leads 
to a natural constraint on the real parts of the Dirac masses Bi{p) and B^{p), 

B,.,{-p) = Bl,{p) . (3-13) 

Using the above properties for Bi and B^ and those for Z\+ and A^ in Eq. ()2-6() . 
we can always restrict the summation over the Matsubara frequencies to a sum over 
only the positive frequencies. 

Substituting the solution of Eq. H3-5|) into Eq. ()3-4() . we obtain the effective 
potential for the vacuum, i.e. at the stationary point. Because the effective potential 
itself is divergent, we subtract the effective potential for the trivial vacuum and 
define 

Vsoi[A+,A-,Bi,Bs] = V[A+,A-,Bi,B3] - F[0, 0,0,0] . (3-14) 

The value of the effective potential in Eq. H3-14() is understood as the difference 
between the energy density of the vacuum at the stationary point and that of the 
trivial vacuum. The true vacuum should be determined by evaluating the value 
of the effective potential: The vacuum with the smallest value of Vsoi is the most 
stable vacuum. The explicit form of the above expression is given in Eq. HA-6|) of 
Appendix El 
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§4. Numerical analysis 

In this section we present the results of the numerical analysis. The only param- 
eter necessary to carry out the numerical analysis is the infrared cutoff parameter 
t/ in the running coupling. The unit of the energy scale in the numerical analysis is 
ylqcd, and it is determined by calculating the pion decay constant /,r for fixed t/ at 
zero temperature and zero chemical potential through the Pagels-Stokar formulap^ 

We use = 88MeV, estimated in the chiral limit as an input. At zero tem- 
perature and zero chemical potential, the dependences of the physical quantities on 
tj have been shown to be small around tf = 0.5 for the running coupling (I) with 
tc = -2.0.1^ For the running coupling (II), we verified that the dependences of the 
physical quantities (mass and condensate) on tf are small around tf = 0.25. There- 
fore, we use tf = 0.5 and tc = —2.0 for the running coupling (I) in Eq. 1)2-13^ and 
tf = 0.25 for the running coupling (II) in Eq. H2-14|) for general T and fj,. From these 
inputs, we obtain ylqcd = 583 MeV for the running coupling (I) and A^^d = 567 MeV 
for the running coupling (II). Later in this section we study the dependence of our 
results on tf. 

We introduce the framework of the numerical analysis in M.ll In ^4.2l we present 
the phase diagram obtained from the SDE in case-1 and the critical exponents of the 
mass gaps and condensates near the phase transition point from the hadronic phase 
to the quark-gluon plasma (QGP) phase and that from the 2SC phase to the QGP 
phase. We investigate the effect of the Debye mass of the gluon in ^4.31 In ^4.41 we 
study the effect of the antiquark contribution by solving the coupled SDEs in the 
three cases introduced in the previous section. 

4.1. Framework of the numerical analysis 

In this subsection we summarize the framework of the numerical analysis. First, 
as discussed below Eq. (|3-13p . we restrict the Matsubara sum to a sum over the 
positive frequency modes using the properties in Eqs. H2-6|) and H3-13() . Second, to 
solve the SDEs numerically, we transform the variable p = \p\ into the new variable 
X. For this transformation, we use a /i- independent transformation in the region of 
small chemical potential (p < /ig) and a //-dependent transformation in the region 
of medium chemical potential {fi > fio). In the region of small chemical potential, 
where the chiral condensate is formed, the characteristic scale of the system is ^Iqcd- 
Here, the dynamical information comes mainly from the region in which p < ylqcd- 
In the high density region, where the diquark condensate is formed, the chemical 
potential fi, in addition to ylqcd; gives an important scale. Here, the dynamically 
important region is that in which p ^ ji. Therefore, we make the transformation 
from p to X as 



X = In [p/ylqcd] for fJ. < fio , 



(4-2) 
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X = \n [p/(3/i)] for n> . (4-3) 

Here, we fix /^o/^qcd = 1/3, because the two schemes of the discretizatfon described 
in Eqs. H4-2() and (|4-3() are connected continuously to each other at this point We 
have confirmed that the two schemes are connected smoothly to each other around 
/X = /xq, as in Ref. 15 ). 

Under the above transformations, the integral over p on the interval [0, oo] is 
converted into an integral over X on the interval [—oo, oo]. In the numerical integra- 
tion, we introduce ultraviolet (UV) and infrared (IR) cutoffs for X, restricting the 
region as X S [Xir, Xjjv]- We discretize the interval of X by introducing Nx evenly 
spaced points: 

X[J] = Xjn + AX-J, J = 0,1,--- ,Nx-l, (4-4) 

where 

AX = h£y^^. (4.5) 

The integration over p is thus replaced by the following summation: 

I dp^ AXY^e-^^'^^- (4-6) 

In the present analysis, for the UV and IR cutoffs, we use 

X: [XirAuv] = [-3.5,2.5] . (4-7) 

For the Matsubara frequency, p^ = (2n + l)vrT, we truncate the infinite sum 
and perform the summation over the 2no modes labeled by integer n in the region 
n S [—no, no — 1]. Throughout this paper, we consider no = 40 and Nx = 40. We 
have checked that these values are sufficiently large for the present purpose. To 
obtain {'tp'ip)iGeV and {Tpip)iGeV, we use Eq. (|2-19|) and Eq. (|2-2Ujl with 

yl = ylqcdexp(Ai7y) for fJ. < flQ , (4-8) 

yl = 3/i exp(A(7y) for ^ > Ho ■ (4-9) 

We solved the SDEs in case-1 with an iteration method. Starting from a set of 
trial functions, we updated the mass functions with the SDE as 



{Si,3oid , ^oid} I Right-hand sides of SDEs (3.6)^(3.9) | {5i,3new , A^^^} . 

(4-10) 



Then, we stopped the iteration when the convergence condition 

"0-1 

72=— no " 



d-^p I \( 5V \V 6V 



(27r)3 4 

no-l „yi 



5[Sf{p)\J \5[Sf{p)\ 

T 22 7^i2|5iold(p)-Slncw(p)|^ + |S3old(p)-S3new(p)|^ 

+2\A+M - Z^+w(p)l' + 2|^-id(p) - ^ncw(p)l'} (4-11 
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was satisfied with suitably small e. In the present analysis, we set e = 10 
case-2 and case-3, we adopted the same method as that used in case-1. 



In 



4.2. Phase structure 

Let us first study the phase strucuture by solving the SDEs (|3-6() - (|3-9|) in case- 
1. The resultant phase diagram obtained by using a running coupling of type (I) 
with tf = 0.5 and Mq 7^ is shown in Fig. ^ In the present analysis, we consider 



2.5 



1.5 



0.5 - 



QGP 



tricritical point 



chiral 
transition 



triple point 



Hadron 



"1 st-order" + 

"2nd-order" x 

"2nd-order-2SC" * 

"tricritical" □ 

"triple" ■ 



diquark 
transition 

2SC 



1 1 1 1 \ — I ' 1 

1 2 3 4 5 

iLi Ik 

Fig. 1. Phase diagram for < T//,r < 2.5 and < fi/ f-n < 5.5. The symbol x denotes the second 
order phase transtion between the hadron phase and the QGP phase, + the first order phase 
transition between the hadron phase and the QGP phase or the 2SC phase, and ^ the second 
order phase transition between the 2SC phase and the QGP phase. Tc=147 MeV, /ic=325 MeV, 
and (T, /i)=(146, 20) MeV at the tricritical point, where the second order phase transition 
changes to first order. At the "triple point", (T, /i) = (110, 192) MeV. Here, the hadron phase 
((V'V) / 0), 2SC phase ((?/>■(/)> / 0) and QGP phase ((?/>^> = (?/>■(/)> = 0) meet. 

the possibility of three phases: the hadron (xSB) phase, the 2SC phase and the 
QGP phase. We determine the true vacuum by evaluating the value of the effective 
potential at the solution of the SDE. The running coupling of type (I) in Eq. (|2-13j) 
used here becomes very large in the infrared region. It exceeds the critical value 
7r/2,*) above which the SDE for the Majonara mass with zero Dirac mass {B = 0) 
provides a non-trivial solution even at T = /i = 0. Then, as we have shown in the 
case n > and T = in Ref. I15|) . even for non-zero but small temperature, the 
2SC vacuum always exits, and in it, a non-trivial Majonara mass is dynamically 

*^ At T — fi — 0, the SDE for A — A" = takes the same form as that for the Dirac mass 
-B, except that an extra factor of 1/2 appears in front of the integration kernel. Because the critical 
value of the coupling in the SDE for B is 7r/4P the extra factor of 1/2 implies that the critical 
value of the coupling in the SDE for A is 7r/2. 
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generated, and the 2SC vacuum is more stable than the trivial vacuum ( the QGP 
vacuum in the present analysis). In the low chemical potential region, the xSB 
vacuum also exists, and in it, the Dirac mass B is dynamically generated, and it is 
most stable among the three vacua. Hence, the true vacuum in the low temperature 
and low chemical potential region is the xSB vacuum. This is natural, because 
the strength of the attractive force mediated by one-gluon exchange between two 
quarks in the color anti-triplet channel is weaker than that between a quark and 
an antiquark in the color singlet channel at T = ^ = 0. When the temperature 
is increased at ^ = 0, the value of the effective potential at the 2SC vacuum goes 
smoothly to zero around T ~ 110 MeV {T/ f^^ ~ 1.3), while the xSB vacuum is still 
the most stable. Up to this region, the value of the chiral condensate (V'V') (per 
flavor) does not change significantly, as we show in Fig. |2 As the temperature is 




Fig. 2. Temperature dependences of the chiral condensate and the Dirac mass gap at /.i = 0. Here, 
the critical temperature of the chiral phase transition from the hadron phase to the QGP phase 
at = is Tc = 147 MeV. 

increased further, the value of the chiral condensate starts to decrease, and it, as 
well as the mass function in the zero-momentum limit (shown by B\l f^^ in Fig. [21), 
finally goes to zero at T/Z^r = 1.67; i.e., 

Tc = 147MeV. (^ = 0) (4-12) 

At the same time, the value of the effective potential in the xSB vacuum also 
smoothly goes to zero, and thus the trivial vacuum appears. Then, at /i = a 
phase transition occurs from the hadron phase to the QGP phase (the chiral phase 
transition), and it is of second order. This is consistent with the results of previous 
analyses performed using the SDE (see, e.g., Refs . IT^ .1^ .ITT I) .[T ^ \ . 
We obtain the tricritical point at 



{Tuiit) = (146, 20) MeV , 



(4-13) 
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(indicated by □ in Fig. For fi > fit the chiral phase transition is of first order 
(indicated by + in Fig. , while for fj, < fit it is of second order (indicated by x in 
Fig. The value of the chemical potential at the tricritical point obtained in the 
present analysis is smaller than those in several other models, such as the instanton 
model, the NJL mode±^'^''^^''^^' and the random matrix model,!^ as well as in the 
analysis done using the SDEjC^'^ but is consistent with that obtained in Eef. II 1)1 . 
which was also obtained through analysis of the SDE. We find that in the framework 
of the SDE, the large difference in the position of the tricritical point in the T-fj, 
plane is caused by the existence of the imaginary part of the Dirac mass. We will 
discuss this point in ^ 

When we increase the chemical potential at T = 0, on the other hand, the 2SC 
vacuum becomes more stable than the xSB vacuum at the critical chemical potential 
fJ-ffiT = 3.69, i.e., 

^c = 325MeV, (r = 0) (4-14) 

before the xSB vacuum becomes less stable than the trivial vacuum. As a result, 
a phase transition occurs from the hadron phase to the 2SC phase, and it is of 
first order. This structure is same as that obtained in Ref. I15() . where the form of 
the gluon propagator was different from the present one, given in Eq. H2-9|) . and a 
running coupling of type (II) in Eq. ()2-14p was used. Let us increase the temperature 
in the 2SC phase. The temperature dependences of the diquark condensate {^Jip) 
(per flavor) and the Majonara mass gap of the quark A~ for fi/ fn = 4 are shown 
in Fig. This shows that the diquark condensate and the Majonara mass gap are 




Fig. 3. Temperature dependences of the diquark condensate and the Majorana mass gap at fi/f-w = 
4. Here, the critical temperature of the color phase transition from the 2SC phase to the QGP 
phase at /i//^ = 4 is fc = W2 MeV. 



stable against change of the temperature in the low temperature region {T/ f.,^ < 0.6, 
i.e., T < 50 MeV). Around T 50 MeV, the condensate and the mass gap start to 
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decrease, and they go to zero at T = Tc = I.IS/tt = 102 MeV. At the same time, 
the value of the effective potential in the 2SC vacuum smoothly vanishes. Then, 
the phase transition from the 2SC phase to the QGP phase occurs, and this color 
phase transition is of second order. This is consistent with the result obtained in 
models with a 4-Fermi interaction.-^' ' However, as we can see from Fig. ^ the 
critical temperature of the color phase transtion obtained in the present analsysis 
is around 100 MeV for any fi satisfying fi < 500 MeV, and this value is larger than 

the values obtained in models based on the contact 4-Fermi interaction.^^''^''^''^ 
This increase of the critical temperature may be caused by the long range force 
mediated by the magnetic mode of the gluon. We will return to this point in the 
next subsection. 

Let us consider the critical behavior of the condensates and the mass gaps near 
the second order phase transition points along the arrows labeled "chiral transition" 
= 0) and "diquark transtion" {n/fn = 4) shown in Fig. ^ First, we consider the 
second order chiral phase transition at zero chemical potential. In Refs. 1^5 ]) . 1131) . IT Hl 
and 139(1 . it was shown, by using SDEs with different kernels, that the scaling prop- 
erties of the chiral condensate and the Dirac mass gap near the critical temperature 
are consistent with those obtained in the mean field treatment. Furthermore, the 
analyses done in Refs. ^ and 138(1 using the SDE with a similar kernel concluded 
that the scaling is consistent with the mean field one, although the numerical errors 
are not clear. Thus, in the present analysis, we assume mean field scaling of the 
condensate and the mass gap, and fit the values of the prefactor and the critical 
temperature*) by minimizing 



where the index i labels the temperature, y is the chiral condensate per flavor —{ijj'il^) 
or the Dirac mass gap in the infrared limit Bi{pQ = p = 0), and F(T, k) is the fitting 
function given by 



Performing the fit for 0.9Tc < T < Tc, we obtain the following best fitted values for 
the prefactor and the critical temperature: 



Biipo = P = 0)~al-— , a~ 15.9/, , T, ~ 147 [MeV] . (4-18) 



We plot the fitting function in Eq. ()4-16(l with the above best fitted values in the left 
panel of Fig. [IJ together with the numerical data obtained from the solution of the 
SDE. This figure shows that the critical behavior of the chiral condensate per flavor, 

*-* We could not precisely determine the critical temperature through consideration of the order 
parameters because of the numerical error introduced by the discretization. 




(4-15) 




(4-16) 




Tc ~ 147 [MeV] , 



(4-17) 
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— {ipip), and the Dirac mass gap, Bi{pQ = p = 0), near fi = are consistent with the 
mean field scaling. 



ln(l-T/1J.) 



ln(l-T/T^) 



Fig. 4. Critical behavior of the condensates and the mass gaps. The left panel shows the critical 
behavior of the chiral condensate per flavor —{tpip) and the Dirac mass gap Bi{po — p = 0) near 
the chiral second order phase transition at n/fn = along the arrow labeled "chiral transition" 
in Fig. while the right panel shows that of the diquark condensate per flavor — (i/ji/j) and 
the Majorana mass gap A~ {po — 0,p — fi) near the diquark second order phase transition at 
fi/ /tt =4 along the arrow labeled "diquark transition" in Fig.Q The solid lines in the left panel 
are the fitting functions given in Eqs. 14- 1711 and 14- 1811 . and those in the right panel are the 
fitting functions given in Eqs. iFT^ and iT^ . 



We next consider the critical behavior of the diquark condensate per flavor 
and the Majorana mass gap on the Fermi surface A~{po = 0,p = /i) near the 
color second order phase transition at ////^ = 4 along the arrow labeled "diquark 
transition" in Fig. ^ Using the same fitting method as used above for the chiral 
transition, we perform the fit for O.OTc < T < Tc- The best fitted values of the 
critical temperature and the prefactor are determined as 

a ~ 11.6/^ , Tc ~ 102 [MeV] , (4-19) 

Y) > « ~ 3.9/^ , fc ~ 102 [MeV] . (4-20) 

We plot the resultant fitting functions together with the numerical results in the 
right panel of Fig. 21 This shows that the critical behavior near the diquark phase 
transition point is consistent with the mean field scaling. 

4.3. Effect of the Dehye mass of the gluon 

In this subsection we study the infiuence of the Debye mass of the gluon on the 
phase diagram in case-1. In Fig. \^we show the change of the phase diagram due 
to the Debye mass. This shows that ignoring the Debye mass does not change the 
following qualitative structure of the phase diagram. There is a tricritical point that 
separates the two branches of the critical line for the chiral phase transition from the 
hadron phase to the QGP phase, the critical line for the second order transition and 



-{ijjip) ~ a ( 1 



1/2 



A (j»o = 0,p = /z) ~ a 1 
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that for the first order transition. The phase transition from the hadron phase to the 
2SC phase is of first order, while the color phase transition from the 2SC phase to 
the QGP phase is of second order. However, the effect of the Debye mass changes the 
quantitative structure. In the case of zero Debye mass, the phase diagram without 
the 2SC phase is same as that in Ref. illjl . However, when we include the 2SC phase 
in drawing the phase diagram, the critical chemical potential becomes smaller than 
that in Ref. HH): At T = 0, /Xc = 460 MeV is obtained in Ref HU, while we obtain 
He = 376 MeV. This is further reduced to fic = 325 MeV by including the Debey 
screening mass into the electric mode of the gluon. Similarly, the critical line dividing 
the hadron phase from the 2SC phase is moved to a region in which the chemical 
potential is smaller by about 15% when the Debye mass is included. Furthermore, the 
inclusion of the Debye mass lowers the value of the critical temperature at the phase 
transition from the hadron phase to the QGP phase by about 15%, while it lowers the 
value of the critical temperature at the phase transition from the 2SC phase to the 
QGP phase by 15—30%. Note, however, that the value of the chemical potential at 
the tricritical point is not changed, although the value of the critical temperature at 
this point decreases by about 15%. These quantitative changes imply that the Debye 
mass of the electric mode plays a role to weaken the attractive interaction between 
two quarks as well as that between a quark and an antiquark. This is reasonable, 
because the Debye mass screens the long range force mediated by the electric mode 
of a gluon. Therefore we expect that the critical temperatures of the chiral and color 
phase transition from the hadron phase and the 2SC phase to the QGP phase may 
be further reduced if the contact interaction, i.e., the short-range force induced by 
the instanton, is the dominant one in this intermediate chemical potential region. 

4.4. Effect of the antiquark contribution 

In this subsection we study the effect of the antiquark contribution by solving 
the following three different types of SDEs: (case-1) Coupled SDEs for A~ and Z\+; 
(case-2) the SDE for A' with Z\+ = 0; (case-3) the SDK for A" with A+ = A' . 
In case-1, the quark mass and the antiquark mass are included properly, while the 
antiquark mass is ignored in case-2. The approximation in case-2 is considered to 
be valid in the high density region. The approximation in case-3 is used in many 
analyses carried out using models with the local 4-Fermi interaction. 

We show the phase diagrams of above three different cases in Fig. |21 In this 
figure, the symbol -|- denotes the phase transtion in case-1, x that in case-2 and >^ 
that in case-3. To make the differences among the three cases clear, we connect the 
data points by a solid curve in case-1, by a dashed curve in case-2, and by a dotted 
curve in case-3. First of all, we should note that the critical line between the hadron 
phase and the QGP phase for h/ fn < 2 is not at all affected by the change of the 
antiquark contribution, because there exist no 2SC solutions. Therefore, in all three 
cases, we obtained the same tricrical point at (T,fj.) = (146, 20) MeV, where the 
critical line is divided into that for the second order phase transition (/i < 20 MeV) 
and that for the first order phase transition (20 MeV < ^ < 180 MeV). Similarly 
to the effect of Debey mass studied in the previous subsection, the effect of the 
antiquark contribution does not change the following qualitative structures of the 
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Fig. 5. Phase diagram for < T/Ztt < 2.5 and < fJ^/ f-n < 5.5 obtained by setting Md = in 
the SDE. For comparison, the phase diagram shown in Fig. Q obtained with the Debye mass 
included is also plotted. Here Tc=U7 MeV, Atc=325 MeV, T^=169 MeV and ^'^=376 MeV. 



phase diagram: In all three cases, the phase transition from the hadron phase to the 
2SC phase is of first order, while the color phase transition from the 2SC phase to 
the QGP phase is of second order. Furthermore, the approximation represented by 
setting = A~ (case-3) results in almost the same phase diagram as in the case 
of the full analysis (case-1), except that the critical chemical potential at T ~ in 
case-3 is slightly smaller than that in case-1. However, the omission of the antiquark 
contribution (case-2) quantitatively changes the critical line between the hadron 
phase and the 2SC phase, as well as that between the 2SC phase and the QGP 
phase. In Fig. |SJ we indicate the region in the 2SC phase in case-1 (but not in 
case-2) by the shaded areas. This clearly shows that the region of the 2SC phase in 
case-2 is smaller than that in case-1, as well as in case-3. The value of the critical 
temperature at the color phase transtion in case-1 becomes larger by about 10% at 
[ij f-ix = 4 than that in case-2. The value of the critical chemical potential at the 
chiral phase transition in case-1 becomes smaller by about 5% at zero temperature 
than that in case-2. 

Let us study the Majorana mass gap and the diquark condensate at zero tem- 
perature. In the remainder of this subsection, we study them not only in the true 
vacuum (i.e., in the case that the 2SC vacuum is the most stable) but also in the false 
vacuum (i.e., in the case that the 2SC vacuum is less stable than the xSB vacuum) 
in order to see the effect of the antiquark contribution more clearly. 

In Fig. [7|we show the dependence on the chemical potential of the quark Majo- 
rana mass gap on the Fermi surface (left panel) and the ratio of the antiquark 
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Fig. 6. Phase diagrams in three different cases for < T/Ztt < 2.5 and < /J./ fw < 5.5. The symbol 
+ denotes the phase transition in case-1, x that in case-2, and ^ that in case-3. The points 
in case-1, case-2 and case-3 are connected by solid curves, dashed curves and dotted curves, 
respectively. The shaded areas are in the 2SC phase in case-1 but not in case-2. 

mass to the quark mass Z\+/Z\^ at po = and p = fi m case-1 (right panel) at 
zero temperature. We note that a nontrivial solution exists in all cases, even at 
= 0, because the running coupling in the infrared region exceeds the critical value, 
7r/2, as discussed in Subsection 14.21 Furthermore, the ratio A~^/A~ is actually 1 
at /i = 0, which is required by the existence of charge conjugation symmetry. As 
a result, in case-3 is equal to that in case-1 at /x = 0. However, the value of 
A~ in case-2 is about half of that in case-1 and case-3 at ^ = 0. When we increase 
the chemical potential, the values of A~ in case-1 and case-3 decrease, while that in 
case-2 first increases and then decreases. The right panel of Fig. [7| shows that the 
ratio Z\+/Z\^ in case-1 decreases as fi increases and reaches a value of about 0.85. 
The left panel of Fig. Q shows that, for n/f-K > 3.5, the values of A~ in all three 
cases are almost same, although the values of differ greatly. In fact, we have 
Z\^(case-1) : Z\^(case-2) : Z\"'" (case-3) ~ 0.85 : : 1. These results imply that the 
antiquark gives a sizable contribution to the quark Majorana mass in the low chem- 
ical potential region, while it becomes negligible for ^//tt > 3.5. This suggests that, 
to form the Majorana mass in the region where the 2SC vacuum is most stable, the 
effect of the Fermi surface is more important than the effect of the large attractive 
force due to the large running coupling. 

Next, we show the dependence on the chemical potential of the diquark conden- 
sates per flavor in three cases in Fig. |H1 This shows that the value of the diquark 
condensate in case-3 {A~^ = A~) is almost same as that in case-1, which is consistent 
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Fig. 7. Dependence on the chemical potential of Z\~(po = 0,p = ^) in three cases (left panel) and 
^^(po = 0, p = /i)/Z\~(po = 0,p = /i) in case-1 (right panel) for < fi/fn < 5.5 and T = 0. In 
the left panel, the data points indicated by the circles, triangles and squares were obtained in 
case-1, case-2 and case-3, respectively. Open symbols (Q, A and □) represent data points in 
the false vacuum (the case in which the 2SC vacuum is less stable than the xSB vacuum), while 
the filled ones represent data points in the true vacuum (the case in which the 2SC vacuum is 
most stable). 



with the fact that in case-1 the value of A'^ is comparable to that of . However, 
there is a large difference between the value of the diquark condensate in case-2 and 
that obtained in case-1 or case-3. Actually, the ratio of the value of the diquark 
condensate in case-1 or case-3 to that in case-2 is roughly 4 in the region of very 
small chemical potential: ('i/'V')case-i / i'^'^) case-2 — 4 for // ~ 0. This can be under- 
stood as follows. Figure 13 shows that the value of A~ in case-1 is about twice that 
in case-2 for ~ 0, and that A'^ is comparable to A~ in case-1 for /i ~ 0. On the 
other hand, the formula for calculating the diquark condensate given in Eq. ()A-5|) 
shows that the diquark condensate consists of quark and antiquark contributions, 
and that the dominant part comes from the ultraviolet region. Thus, the quark and 
antiquark contributions to the diquark condensate in case-1 are each about twice 
the quark contribution in case-2, and therefore the ratio of the value of the diquark 
condensate in case-1 or case-3 to that in case-2 is about 4 near ^ = 0. Now, how 
about the ratio in the region of intermediate chemical potential {fJ^/ fn ^ 3.5)? If the 
effect of the Fermi surface were negligible, the same argument as above would lead 
us to conclude that the values of diquark condensate in case-3 and case-1 are roughly 
twice that in case-2. However, Fig. |H1 shows that the ratio is actually smaller than 2, 

and it becomes about 4/3 at jj/ /t, = 5: (•i/'^)case-i / i'^i') case-2 — ^/^ ^'^^ t^l f-^ — ^• 
This implies that, as expected, the effect of the Fermi surface plays an important 
role in the region of intermediate chemical potential. 

The results we have found in this subsection to this point (Figs. inHHl) imply that 
the antiquark gives a non-negligible contribution, and that the approximation A~ = 
A^ is sufficiently valid to allow a study of the phase diagram, the quark Majorana 
mass gap and the diquark condensate in the region of intermediate temperature 
and intermediate chemical potential, where the chiral phase transition occurs. The 
same approximation may be sufficiently valid to allow investigations of other physical 
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Fig. 8. Dependence on the chemical potential of the diquark condensates per flavor for < /i//,r < 
5.5 and T — 0. The data points indicated by the circles, triangles and squares were obtained in 
case-1, case-2 and case-3, respectively. The open symbols (Q, A and □) represent data points 
in the false vacuum (the case in which the 2SC vacuum is less stable than the xSB vacuum), 
while the filled ones represent data points in the true vacuum (the case in which the 2SC vacuum 
is most stable). 



quantities such as the number density, as well. 

Let us compare the results in the ratio of the critical temperature to the Ma- 
jorana mass gap at zero temperature with the BCS result, Tc/A ~ 0.567, obtained 
in the high density region.^'EHJ In Fig. |21we show the dependence on the chemical 
potential of the ratios of the critical temperature to the quark Majorana mass gap 
on the Fermi surface at zero temperature. There are only small differences among 
the data points obtained in the three cases. From this figure, we see that this ratio 
approaches the BCS value, 0.567, when we increase the value of the chemical po- 
tential, and that it is already close to the BCS value in the intermediate chemical 
potential region: 

Tc/A^ 0.5 , for ^ ~ 450 MeV . (4-21) 

Finally, we check the dependence of the phase diagram on the model used by 
changing the infrared regularization parameters for two types of running couplings 
given in Eqs. H2-13() and (|2-14() . We show the phase diagrams for the running coupling 
of type (I) with tj = 0.4, 0.5 and 0.6 in the left panel of Fig. and the phase 
diagrams for the running coupling of type (II) with tj = 0.20, 0.25 and 0.30 in the 
right panel of Fig. 1101 These figures show that the critical line for the phase transition 
from the hadron phase to the 2SC phase, as well as that from the 2SC phase to the 
QGP phase, has a weak dependence on the infrared regularization parameter tf. For 
a larger value of tf, the phase transition from the hadron phase to the 2SC phase 
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Fig. 9. Dependence on the chemical potential of the ratios of Tc to Z\~(p4 = 0-p= fj,)\^_^ for 
< n/fn < 5.5. The data points indicated by the circles, triangles and squares were obtained 
in case-1, case-2 and case-3, respectively. The open symbols (O; ^ a-nd □) represent data points 
in the false vacuum (the case in which the 2SC vacuum is less stable than the xSB vacuum), 
while the filled ones represent data points in the true vacuum (the case in which the 2SC vacuum 
is most stable). 

occurs at a smaller chemical potential, and that from the 2SC phase to the QGP 
phase occurs at a smaller temperature. Contrastingly, the critical line for the chiral 
phase transition from the hadron phase to the QGP phase is very stable against 
changes in for both types of the running coupling. Furthermore, the running 
coupling of type (II) with tf = 0.25 gives almost the same critical line for the chiral 
phase transition as that of type (I) with tj = 0.5. These results imply that the phase 
structure in the region of small chemical potential is very insensitive to the infrared 
regularization of the running coupling. 

§5. Summary and discussion 

We studied the phase structure of hot and/or dense QCD by solving the Schwinger- 
Dyson equations for the Dirac and Majorana masses of the quark propagator with 
the improved ladder approximation in the Landau gauge. We found that there exists 
a tricritical point at (T , /u) = (146, 20) MeV, which divides the critical line between 
the hadron phase and the QGP phase into the line of a second order phase transition 
[between (T , fi) = (147, 0) MeV and (146, 20) MeV] and that of a first order phase 
transition (20 MeV < fi). Our result implies that at the second order phase transi- 
tion at /X = 0, the scaling properties of the Dirac mass and the chiral condensate are 
consistent with the mean field scalings. 
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Fig. 10. Phase diagram for several choices of the infrared cutoff parameter tf. The left panel 
plots the critical lines for tj = 0.4, 0.5 and 0.6 with the running coupling of type (1), given in 
Eq. 12-131 1. The right panel plots the critical lines for t/ = 0.2, 0.25 and 0.3 with the running 
coupling of type (II), given in Eq. H2-14|l . 



The phase transition from the hadron phase to the 2SC phase was found to be 
of first order at non-zero temprature, as we have obtained at zero temperature in 
previous analysis.^ Furthermore, we found that the color phase transtion from the 
2SC phase to the QGP phase is of second order, with the scahng properties of the 
Majorana mass and the diquark condensate consistent with the mean field scalings. 
The resultant critical temperature of the color phase transition is on the order of 
100 MeV, which is about twice the value obtained in analyses carried out with models 
based on the contact 4-Fermi interaction (see, e.g., Refs. [TH))MI)) J^T |l ] ^ ). We believe 
that this increase of the critical temperature may be caused by the long range force 
mediated by the magnetic mode of the gluon. 

In the present paper we performed analysis with the imaginary part of the Dirac 
mass included, because in the SDE at non-zero chemical potential, the imaginary 
part, which is momentum dependent (an odd function in po), is inevitably generated 
in the hadron {xSB) phase. However, some analyses of the SDE do not include 
the imaginary part, and analyses done using the local 4-Fermi interaction model do 
not generally include the imaginary part, because a leading order approximation is 
used. We found that the most noteworthy feature of the analysis that includes the 
imaginary part of the Dirac mass regards the position of the tricritical point in the 
T-// plane. When we use the SDE without the imaginary part of the Dirac mass 
(i.e., Im[S(p)] = 0), the tricritical point is at (r,/i) = (124,210)MeV [and the triple 
point is at (T, //) = (111, 243)MeV]. The value of /i at this tricritical point is close to 
the values [// ~ O(100)MeV] obtained from analyses carried out using models with 
the local 4-Fermi interactioiP*^^''^^ and those carried out using the SDEs 
with a momentum dependence of the mass function assumed.^^'' The result here 
implies that, in the SDE analysis, including the imaginary part causes the tricritical 
point to move to a position of smaller chemical potential: (T, ^) = (124, 210)MeV 
{T,fi) = (146, 20)MeV. Therefore, we believe that the imaginary part of the Dirac 
mass is important and should be included in the SDE at finite chemical potential. 

We studied the effect of the Debye mass of the gluon, and showed how the phase 
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structure is changed. When we ignore the Debye mass, the critical temperature for 
the color phase transition from the 2SC phase to the QGP phase is increased by about 
15—30%. In addition, the critical temperature Tc for the chiral phase transition from 
the hadron phase to the QGP phase and the critical chemical potential for the 
phase transition from the hadron phase to the 2SC phase are also increased by about 
15%. 

We examined the effect of antiquark contribution. Our results show that the 
antiquark Majorana mass gap, Z\^, is comparable to the quark one, Z\^, in all 
the chemical potential regions that we studied: 1 > /A~ > 0.85 for < n < 
500 MeV. As a result, the approximation represented by ignoring the antiquark mass, 
which is generally considered to be good at very high density, results in quantitative 
differences in the phase diagram, the value of the Majorana mass gap, and that of 
the diquark condensate. 

By contrast, analysis employing the approximation represented by setting = 
A~ , which is often adopted in analysis using models with the local 4- Fermi inter- 
action, gives almost the same results for the phase diagram, the quark Majorana 
mass gap and the diquark condensate as those obtained using the full analysis, in 
spite of the fact that the value of is about 15% smaller than that of A~ in the 
full analysis. Presumably, in regions of intermediate temperature and intermediate 
chemical potential, it is also sufficient to set the antiquark Majorana mass equal to 
the quark one {A^ = A^) for the purpose of studying other physical quantities, such 
as the number density. 

As in previous analysis done at T = 0,^^ we sought the mixed phase, in which 
both the chiral condensate and the diquark condensate exist, by solving the coupled 
SDEs for Majorana and Dirac masses, starting from several initial trial functions. 
However, we could not find such a solution in the present analysis. 

It is important to investigate the QCD phase structure by studying phenomena 
in compact stars and heavy ion collisions. In this paper, we find the color phase 
transition to be of second order. If this is the case, there is a stronger possibility for 
the existence of a pseudogap phase, as discussed in Ref . 0T|) . 

In this paper we ignored the Landau damping of the magnetic mode of the 
gluon, because the approximated form adopted in the previous analysi^^ may not 
be valid in the region of low chemical potential. It may be interesting to study 
the effect of Landau damping by including it as a hard thermal and/or dense loop 
correction without further approximation, although we may have to perform the 
angular integration of the SDE numerically in such a case. 

We assumed that the Landau gauge causes no deviation of the wave function 
renormalization of the quark propagator from 1, even at non-zero temperature and/or 
non-zero chemical potential. In order to keep the wave function renormalization 
equal to 1, we may have to introduce a non-local gauge fixing term, as in Ref. I28|) . 
at T = ^ = 0. We expect that this will not change the qualitative structure of 
the present results. In this paper, we have the QCD scale Aqcd ~ 600 MeV, which 
is larger than the value determined from the experimental value of in the high 
energy region, although the value of the pion decay constant is set to a value 
consistent with experiments. In Ref. 14 Ij) . an effective coupling that includes higher 
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order corrections is used to show that the values of A^cd and /tt become consistent 
with experiments. It would be interesting to use such a running coupling to analyze 
hot and/or dense QCD. 
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Appendix A 

Condensates, Effective Potential and Schwinger- Dyson Equation 



In this appendix we give explicit forms of condensates, the effective potential 
and SDEs. The integration kernels in the SDEs are different from those in Ref. I15|l . 
because the forms of the gluon propagators used in the SDEs are different. 

The explicit forms of the chiral condensate and the diquark condensate are given 

by 



Nr 



F{p,Bi,A) 



d^p 
(2^ 

{po-tif-p^ - {Bi{-p)V ]Bi{p) - A+{p)A-{p)B^{-p) 



1 



F(p,53,^ = 0) 

I.T^a 



= 4(iV,-l)iV^T^ 



[p^-l^f -p'' -{Bz{-p)Y]B:,{p) 



d^p 



(A-1) 



(27r)3F(p,Si,Z\)2 



{po? - iP + l^f - {A+{p)}' - \B^ip)\^ }A-ip) 



+ (Po? - iP-f^f - {^-(P)}' - \Bi{p)\^ \A+{p) 



(A-2) 



where F is defined as 



F{p, B, A) = [{po + i^f -p^- {Bip)}'] [{po - fif -p'- {B{-p)}^] 

-[{pof - ip-^rwA+ipr - [{po? - {p+finA-ip)\' 

+\A+ip)\^\A-{p)\^ + 2Bip)Bi-p)A+ip)A-ip) . (A-3) 

The chiral condensate in Eq. ()A-1|) takes the following form when we set A = and 
Bi = Bs = B: 
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AN^NfTY, f 



B{p) 



(27r)3(po + //)2-p2_{iJ(p)}2 



(A-4) 



where Nc = 3 and Nf = 2. The diquark condensate in Eq. HA-2|) takes the following 
form when we set Bi = i?3 = 0: 



4iN,-l)NfTY^ J" 



(2^ 



A-{p) 



+ 



A+{p) 



(po)2-(p-m)2-{/^-(p)}2 ipoy-ip + fir-{A+ip)y 



(A-5) 



In this expression, the first term in the square brackets is the contribution from the 
quark, and the second term is that from the antiquark. 

The explicit form of the effective potential H3-14() is given by 

= l'[zl+,/\-,Bi,B3]-F|0, 0,0,0] 



(fp 

(27r)3 



2 In 



F{p,Bi,A) 



+ ln 



dPp 



[{po + - p^] [(po - ^J'Y - 

F{p,B^,Q) 

[(po + - p^] [{po - i^Y -p^ 

2 



[(Po - - r - {B,{-p)Y][{p, + isY - p'] 



+ 



F{p,Bi,A) 

+ [iPo + -P'- {Bi{p)}'][{po - - P'] 
-[{Pof - ip + f^nA~\' - [{pof - ip-f^fWA+l' 

[{po - m)' -p'- {B3{-p)}'] [(po + f^f - p'] 



Fip,B3,0) 



+ [ipo + l^f -P'- {B3{p)}'][ipo - l^f - P'] 



(A-6) 

Substituting the expression for Spii into the SDEs for Bi and -B3 in Eqs. (|.S-6j) 
and p-7|) . we obtain 



no-l 

B,{p) = -T 



dqq^ 



(2^)- 



■27ra= 



n=— no 

xKo{q4,P4,q,p) 



5 F+{q,Bi,A) 1 

[6 F{q,Bi,A) 2{iq^ + ^iY-(p-{BMY 



(A-7) 
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no-l 
n=— no 



dqq^ 



(27r)3 
xKo{qA,PA,q,p) 



27rao 



F+iq,Bi,A) 1 ^3(^^ 

F{q,Bi,A) +3(% + ^)2-g2-{S3(9)}2 



(A-8) 



where 



F+{q,B^,A) = B^{q)[{iq^ - iif - f - {B^{-q)f] - B^{-q)A+{q)A-{q) , 



with 54 



(A-9) 

(2n + 1)itT and p4 = — ipo- The integration kernel Kq is given by 



Ko{q4,PA,q,p) =-i / dnDf,^{q - p)tT{j>'A^j'') 



47r, |g + p|2+a;2 27r ^ |g + p|2 + a;2 + 2M?, 

— log -rz =12- 2 + ~ TZ Zfn- 2 I o)i#2 



where 

^ = k4 -P4| • 

When we set A~ = Z\+ = 0, we have 

F+{q,Bi,A = 0) Bi{q) 



F{q,B,,A = 0) (ig4 + ^)2_g2_{5^(^)}2 



(A-IO) 
(A-ll) 

(A-12) 



If we further set Bi = B^, the two equations in Eqs. (|A-7|) and ()A-8|) become identical. 
This implies that Bi = B^ is actually a solution of the SDEs for A^ = Z\+ = 0. 

Next, substituting the expression for Sfi2 into the SDEs for A^ and Z\+ in 
Eqs. and (|3^ . we obtain 



no-l 



n=— no 



(2^)- 



■27rao 



■ _ 2 _ 2 G.{q,B,,A) 

Ki{q4,PA,q,P) ■ - ■ -pTr — 5 — + K2{q4,P4,q,P) ■ - 



3 F(g,5i,zA) 



3 F{q,Bi,A) 

(A-13) 



Z\+(p) 



no-l 

^ E 

n=— no 



■27rao 



Kiiq4,P4,q,p) ■ - 



2 



Fiq,B,,A) 



+ K2{q4,P4,q,p) ■ - 



2 G+{q,Br,A) - 

3 ■ F((?,Bi,zA) J ' 
(A-14) 
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where 

B^,A)= BMBi{-q)A+{q) + [{q^f + {q + + |Z\+(<?)|2]Z^-(g) , 

(A-15) 

A) = B,{q)Br{-q)A-{q) + [{q^f + (g - + \A- {q)\^]A+ (q) . 

(A-16) 

The integration kernels Ki and are given by 



Ki{qA,P4,q,p) 



dOD^^iq-pMA^j^^AjY) 



+ 



2q^p'^Lo'^ {q — p)"^ 2q^p^ 



In + ^_o_n S ^ h 4gp - Hn ■ 



(g — py + a;^ 



^ {iq+p)^ + 2Ml + (^^}{(2M|) + LO^ + {q - pfuo^} (q + pf + 2MI + u 



2M2 + cj2 (g - p)2 



2M2 (2M|, + 



{(g + p)^ + ^^}{{q - P? + + 



2M2 



(g - p)2 + 2M2 + 

(A-17) 



K2{q4,P4,q,p) = -i j dQDp^uiq - p)tMpl^ A^-f""^ 



TT{q^-p^f {q + pf 



TT \{q^ - p^) 



;2\2 



[q — p)2 + u;2 



2g2p2 

^ {(g - + 2M2 + a;^}{(2M2 + ^^2)2 ^ ^ j^)2^2| + ^)2 ^ + a;2 



2M2, + u;2 {q- pf 



2M2 (2M2, + 



+ pf + '^^lil^ - P? + (g + p)^ + w^"' 



2M2, 

We note that for i?i = 0, we have 
G±(g,Si =0,Z\) 



{q — pY + 

-A^{q) 



F{q,B,=0,A) (g4)2 + (^-^^)2 + |^T(g)|2 



(g - p)2 + 2M2 + (j2 

(A.18) 



(A.19) 



Thus, the SDEs ()A-13|) and ()A-14|) are identical to the well-known forms given 
in, e.g., Ref.lHl). 



References 



1) G. E. Brown and M. Rho, Phys. Rep. 363 (2002), 85, |hep-ph/0103102j Phys. Rep. 269 

(1996), 333, hep-ph/9504250 
F. Wilczek, hep-ph/0003183 
R. D. Pisarski, hep-p h/ 9503330] 



S. Takagi 



2) S. P. Klevansky,"hep-ph/9810399" 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994), 221, |hep-ph/9401310| 

3) K. Rajagopal and F. Wilczek, hep-ph/0011333 

M. G. Alford, Ann. Rev. Nucl. Part. Sci. 51 (2001), 131, |hep-ph/0102047| 

4) K. Kanaya, hep-ph/0209116j 
Z. Fodor, hep-lat/020910lp 

5) V. A. Miransky, Dynamical symmetry breaking in quantum field theories (Singapore, Sin- 
gapore, World Scientific, (1993), p. 533). 

T. Kugo, "Basic concepts in dynamical symmetry breaking and bound state problems," 
KUNS-1086 Lecture delivered at 1991 Nagoya Spring School on Dynamical Symmetry 
Breaking, Nakatsugawa, Japan, Apr 23-27, 1991. 

C. D. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. 45 (2000), SI, nucl-th/0005064 

D. T. Son, Phys. Rev. D 59 (1999), 094019, hep-ph/9812287 

D. K. Hong, V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Rev. D 61 
(2000), 056001, [ibid. D 62 (2000), 059903, Erratum], hep-ph/9906478 
T. Schafer and F. Wilczek, Phys. Rev. D 60 (1999), 114033, hep-ph/9906512;| 
Y. Taniguchi and Y. Yoshida, Phys. Rev. D 55 (1997), 2283, hep-ph/9512375"^ 
D. Blaschke, C. D. Roberts and S. M. Schmidt, Phys. Lett. B 425 (1998), 232, 
nucl-th/9706070 

M. Harada and A. Shibata, Phys. Rev. D 59 (1999), 014010, hep-ph/9807408l 
A. Barducci, R. Casalbuoni, S. De Curtis, R. Gatto and G. Pettini, Phys. Rev. D 41 
(1990), 1610. 

A. Barducci, R. Casalbuoni, G. Pettini and R. Gatto, Phys. Rev. D 49 (1994), 426. 
O. Kiriyama, M. Maruyama and F. Takagi, Phys. Rev. D 62 (2000), 105008, 
hep-ph/0001108 

O. Kiriyama, M. Maruyama and F. Takagi, Phys. Rev. D 63 (2001), 116009, 
hep-ph/OlOlUO 

M. Harada and S. Takagi, Prog. Theor. Phys. 107 (2002), 561, hep-ph/0108173 
K. Rajagopal and E. Shuster, Phys. Rev. D 62 (2000), 085007, hep-ph/0004074 
R. D. Pisarski and D. H. Rischke, Phys. Rev. D 61 (2000), 051501, nucl-th/9907041 Phys. 
Rev. D 61 (2000), 074017, nucl-th/9910056 

J. Berges and K. Rajagopal, Nucl. Phys. B 538 (1999), 215, |hep-ph/ 9804233 
R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Ann. of Phys. 280 (2000), 35, 

hep-ph/9904353 ^ 

H. Mishra and J. C. Parikh, Nucl. Phys. A 679 (2001), 597, hep-ph/0003019 
T. M. Schwarz, S. P. Klevansky and G. Papp, Phys. Rev. C 60 (1999), 055205, 
nucl-th/9903048 

M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, Phys. Rev. D 65 (2002), 091504, 
nucl-th/0111022 

M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, hep-ph/0207255j 

B. Vanderheyden and A. D. Jackson, Phys. Rev. D 62 (2000), 0940 fn; |hep-ph/0003150l 
K. i. Kondo and K. Yoshida, Int. J. Mod. Phys. A 10 (1995) 199, hep-th/9304018| 
A. Niegawa, Mod. Phys. Lett. A 17 (2002), 303, hep-th/0203045 

M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B 422 (1998), 247, hep-ph/9711395 
R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81 (1998), 53, 
hep-ph/9711396 

M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537 (1999), 443, hep-ph/980'4403l 
T. D. Fugleberg, hep-ph/0206033 

T. Kugo and M. G. Mitchard, Phys. Lett. B 282 (1992), 162. 

M. Le Bellac, "Thermal Field Theory" (Cambridge University Press, Cambridge, England, 
1996). 

K. Higashijima, Phys. Rev. D 29 (1984), 1228. 
V. A. Miransky, Nuovo Cim. A 90 (1985), 149. 

K. I. Aoki, M. Bando, T. Kugo, M. G. Mitchard and H. Nakatani, Prog. Theor. Phys. 84 
(1990), 683. 

K. I. Aoki, T. Kugo and M. G. Mitchard, Phys. Lett. B 266 (1991), 467. 
J. M. Cornwall, R. Jackiw and E. Tomboulis, Phys. Rev. D 10 (1974), 2428. 

C. De Dminicis and P.C. Martin, J. Math. Phys. 5 (1964), 31. 



Phase Structure of Hot and/or Dense QCD 



31 



34) H. Abuki, T. Hatsuda and K. Itakura, Phys. Rev. D 65 (2002), 074014, hep-ph/0109013" 

35) J. C. Bloch, C. D. Roberts and S. M. Schmidt, Phys. Rev. C 60 (1999), 065208, 
|nucl-th/9907086 

V. A. Miransky, I. A. Shovkovy and L. C. Wijewardhana, Phys. Lett. B 468 (1999), 270, 
|hep-ph/9908212 

N. Evans, J. Hormuzdiar, S. D. Hsu and M. Schwetz, Nucl. Phys. B 581 (2000), 391, 
|hep-ph/9910313 

lA. Matsuzaki, Phys. Rev. D 62 (2000), 017501, hep-ph/991054l] 

36) H. Pagels and S. Stokar, Phys. Rev. D 20 (1979), 2947. 

37) J. Gasser and H. Leutwyler, Ann. of Phys. 158 (1984), 142. 

38) T. Ikeda, Prog. Theor. Phys. 107 (2002), 403, hep-ph/0107105l 

39) A. HoU, P. Maris and C. D. Roberts, Phys. Rev. C 59 (1999), 1751, "nucl-th/9808076' 

40) D. Bailin and A. Love, Phys. Rep. 107 (1984), 325. 

41) M. Hashimoto and M. Tanabashi, hep-ph/0210115 



